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Abstract. We derive a time-convolutionless master equation for the spin-boson 
model in the weak coupling limit. The temporarily negative decay rates in the 
master equation indicate short time memory effects in the dynamics which is explicitly 
revealed when the dynamics is studied using the non-Markovian jump description. The 
approach gives new insight into the memory effects influencing the spin dynamics and 
demonstrates, how for the spin-boson model the the co-operative action of different 
channels complicates the detection of memory effects in the dynamics. 
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1. Introduction 

The spin-boson model describes a two-level system interacting with a bosonic 
environment. It is widely studied, mostly because of its many applications in the field 
of quantum information theory and chemistry. The model can describe for example 
a simple two- well potential [1], double quantum dots [2j and different biomolecular 
systems [31 Hj. Because the spin-boson model is not exactly solvable many approximative 
approaches have been used in order to unravel the dynamics of the spin-boson system. 
The model has been studied using perturbation theory in the system-environment 
coupling with the Markov approximation [5], [6] and projection operator methods [7J. In 
reference [Ij an approximative scheme called the interacting blip approximation (NIBA) 
has been introduced. It is a perturbative approach in the tunneling matrix element 
based on path integral methods. There also exists an exact approach to the problem 
using Feynman- Vernon path integral method [HI [9]. However, this presents only a formal 
solution. 

Understanding the short-time dynamics of the spin-boson model is crucial in 
understanding quantum computing systems [TOl [TT] . The Markovian approximation 
is valid only for long times and thus the non-Markovian description of the dynamics is 
needed. NIBA fails for large tunneling matrix element and cannot therefore be applied 
to biomolecular systems [S, Hj. Many projection operator methods and path integral 
methods are difficult to treat and do not necessarily give much physical insight to the 
problem. The main motivation of this paper is to write down a time-convolutionless 
master equation describing the dynamics of the spin-boson model and to study the 
equation via non-Markovian quantum jumps. This description gives a clear picture of 
the processes involved in the dynamics and reveals the role of memory effects in the 
evolution of the spin. 

In section [2] we derive a time-convolutionless master equation describing the 
reduced dynamics of the spin-boson model by using the time-convolutionless projection 
operator method fl2]. The master equation is derived in second order and the secular 
approximation is performed. The derived master equation is then studied in section [3] 
via the non-Markovian quantum jump (NMQJ) approach [T3l[Tll[T5] and the conclusions 
are presented in section HI 

2. Derivation of the master equation 

The Hamiltonian describing the spin-boson model is 



Here, = |e) (e| — \g) {g\, = \e) {g\ + \g) (e| in the eigenbasis of the two-state 
system {|e) , \g)}, e is the energy bias between the two states, A is tunneling amplitude 
between the states, and Un and a„ are the energies and annihilation operators of the 
corresponding modes of the bosonic environment respectively. The tunneling amplitude 
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A is assumed to be a real number. The parameter A„ describes the couphng strength 
between the two-state system and the nth mode of the reservoir. The couphng of 
the system to the environment is a bihnear couphng through operator cr^ ,i.e., the 
environment is sensitive to the value of the operator a^. 

We start the derivation of the master equation by writing down a general form of 
the second order time-convolutionless master equation [12] 



dt' 



ps{t) = - / tie[Hi{t), [Hi{t'),ps{t)^ps]]dt'. 



(2) 



Here, the notation A{t) stands for an operator A written in the interaction picture. Now, 
if one writes a decomposition for the interaction Hamiltonian of the form Hj = S ® E 
in terms of the system eigenoperators S^i, i.e. Hi{t) = J^ui ^"^'^^ E(t), and performs 
the secular approximation, ([2]) may be written in the form 

-psit) = -i[HLsit),psit)]+Vit)ipsit)), (3) 



with 



and 



miPsit)) = E7c.(t) \s^ps{t)Sl - \ {Sls^rpsit)} 



(4) 



The time-dependent coefficients 7a; (i) and X^iit) in 
expression 

r.(t) = 



(5) 

and ([5]) are determined from the 

<lt'e'^''ii£[E\t)E{t~t')ps} 

<lt'e'''''iie[E\t')E{Q)pe} (6) 

by the relation Ti^{t) = \juj(t) + i\^(t). 

The system eigenoperators for the spin-boson model written in the system 
eigenbasis {ip+j'ip-} are 

^0 = t; — Si^g = —- — (T_, 5'_a;o = ~n — '^+' (''') 



2uJo 



2uJo 2ujq 
a_ = (tp^l, a+ = and ujq = Ve^ + 



where az 

is the system eigenenergy. In the limit of continuous spectrum of the environment and 
for zero temperature one can write the coefficients 7a; (i) in dl]) as 



dt' diu'J{iu')cos[{iu-iu')t'], 



where J{uj) is the spectral density. For an Ohmic spectral density J{oj) 
[T2] . where Uc is the cutoff frequency, the coefficients are 



7 cue" 



(8) 

-Lj/uJc 



[uJct cos{ijjQt) =F sin(co'ot)] 



-|- auJoe 



[Si(z)] T $5 [Ci(z)] ± 



TT 



7o(t) 



(9) 
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where Ci{z) is the cosine integral, Si{z) is the sine integral and z = Uot + i^. The decay 
rates 7<^o(^) and 7_(^(,(t) have oscillatory behavior and they get temporarily negative 
values resembling memory effects in the reduced dynamics, while the decay rate 70 (t) 
remains always positive. For simplicity we ignore the Lamb shift and restrict to the 
regime Uc <^ ujq, where the secular approximation is valid in the non-Markovian regime. 

The solution to the master equation presented in ([3]), ©and ([9]) can be written in 
terms of the density matrix elements as 



p{t) = M(t)p(0), 



(10) 



with 



M{t) 



Here, 



git) 
fit) 

vit) 

m 
m 



( git) 

e~^W 



\l-git) 



fit) + e-^W 
A2 






-C(i) 



fit) \ 




1 - fit) ) 



7-.o(^)e-''^*Mt, 



A 



i^loit) + ^ilu^M + 7--„ W))dt. 



(12) 



The matrix in equation (|TOl) acts on p(0) as on a column vector consisting of the density 
matrix elements. 



3. Jump description of the dynamics 

The master equation derived in the previous section can be studied with the NMQJ- 
method [13] , which is a generalization of the Monte Carlo wave function method [16] to 
the non-Markovian regime. The key increment in this description is that whenever the 
coefficients 7a; (t) get negative values the direction of the jump process is reversed for 
each wave function realization in the corresponding channel. The density matrix for an 
ensemble of wave functions with N members at time t is 

psit)=j:^\Mt)){Mt)\, (13) 

where Nait) is the number of ensemble members in state |0a(t)) at time t. When 
the decay rates are positive the jumps taking place in a single wave function history 
can be determined by the Monte Carlo wave function method [16j. When a decay 
rate becomes negative, the reversed jump operator in the corresponding channel j is 
Di_^At) = |0„'(t))(0a(t)|, With |0„(t)) = C,_\<f)'^{t))/\\C,_\<f)'^{t))\\. The reversed 
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jumps can take the wave functions to superposition states destroyed earlier by jumps 
and they therefore resemble the memory effects in the dynamics. The memory effects, 
by recreation of superpositions, allow a reversal of the decoherence process which we call 
here recoherence. The reversed transition occurs during a time step 6t with probability 

iV„,(t) 



At) 



(14) 



Let us study the dynamics of the spin-boson model via the jump description. 
We concentrate here on giving a clear picture of the effect of the reversed jumps on 
the dynamics and will not perform the actual unraveling of the master equation ([3]). 
The intention is to give insight into the origin of memory effects by describing the 
system dynamics in terms of non-Markovian jumps. The plots describing the dynamics 
presented in this section are obtained by numerically evaluating (ITT!) . 

The spin-boson model involves three jump operators Ci = cr_, C2 = o-+ and 
C3 = cr^. The operators Ci and C2 generate jumps that affect the populations of 
the system eigenstates and the operator C3 produces phase flips. The corresponding 
decay rates are 71 (t) = ^l^oit)^ 72 (i) = ^7-<^o(^) and 73 (t) = ^7o(i) (cf- © and 
(Q ) , which are plotted in figures [T|2] and [31 The first two decay rates have temporarily 
negative values while the third decay rate remains always positive. This indicates that 
in channels Ci = o"_ and C2 = 0"+ there occurs reversed jumps resulting into memory 
effects, but that the phase flip channel C3 = cr^ lacks memory. 



ri(t) 
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Figure 1: The decay rate 7i(t) in units of Uc for a = 0.01, ojq/uJc 
e/A = l/(2v^). 
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Figure 2: The decay rate 72 (t) in units of Uc for a 
e/A = l/(2v^). 
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Figure 3: The decay rate 73 (t) in units of Uc for a 
e/A = 1/{2VS). 
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Assume now that the system is initially in a pure state iV'o(O)) = + I'lp-)), 

i.e., all the ensemble members are initially in the same state. We use the notation \iIjq) 
for the state vector exposed to only deterministic evolution. When e 7^ 0, there exists 
a phase flipped counter part {ipQ^) for the state {ipo), into which the single realization 
can jump through channel C3 = o"^. The possible states and jumps are demonstrated 
in figure HI 




Figure 4: The states involved in the effective ensemble f|T3l) . The possible jumps in the 
dynamics are presented as arrows between the states. Here, \ipo) = ao\4'+) + Poli^-) 
and \^P^'')=ao\^P+)-/3o\^p-). 



In figure O the coherence is plotted for a fixed ratio Wq/uJc and varying values of 
e/A. One can see that as e/A varies, the behavior of the coherence changes in such 
way that for small values of e/A there is recoherence but as the ratio gets larger the 
system decoheres with no revivals of coherence. Whether the decay rates 71 (t) and 
72(t) get negative values depends only on the ratio Uq/uc, so for large values of e/A the 
system will decohere although there are reversed jumps in the process. Let us study the 
real part of the coherence Re [p±] = cuiPi, where l^po) = Oii + (3i Jumps 
iV^o) \i)+),\i)o) -> IV"-), IV'™) 1^^+), iV'o") and phase flips reduce the 

number of terms in the sum presenting the coherence and therefore induce decoherence. 
Reversed jumps \ipo) !'?/'+), |V'o) ^ |^-) and IV"™) ^ 1"^+)) iV'o^) ^ IV"-) produce 
new terms into the sum, but if there is an equal amount of reversed jumps to iV'o) and 
I'?/'™) the new terms in the sum cancel each other and no recoherence takes place. So 
the appearance of recoherence does not only depend on the number of reversed jumps. 
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Figure 5: The coherence for ^ = 10 and = (black hne), ^ = 10 and = 0.255 
(dark gray hne), ^ = 10 and = (hght gray hne). Here, the couphng constant is 
a = 0.01 and the initial state is \ipo{0)) = "^(1"^+) + l"^-))- 

but also on the proportion of the number of jumps into different target states. 

To find out for which values there is recoherence let us study the quantity 
No(t) — N^^it). Here, No(t) is the number of ensemble members in state \ipo) and 
A^™(t) the number of elements in state The real part of the coherence is an 

mcreasmg function, when No{t) - N^^{t) is increasing, and a decreasing function when 
No{t) - N^^{t) is decreasing. The regions, where there is recoherence for a fixed value 
of ujq/oJc can be seen in figure El One can see that as the energy bias e becomes larger 
compared with the tunneling amplitude A the probability of a phase flip increases. As 
more phase flips occur a reversed jump to IV"™) gets as probable as a reversed jump to 
state which then causes the system to decohere. 
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Figure 6: The time regions for recoherence for different values of In light brown 
areas there is recoherence. Here, ^ = 10, a = 0.01 and \iPq{Q)) = ^(|'?/'+) + The 
values of used in figure [5] are marked in the plot. 

The study above illustrates how for the spin-boson model there are cases where there 
are no revivals of coherence, but there occurs reversed jumps indicating memory effects 
on the level of single wave function realizations. There are no revivals of coherences due 
to the co-operative action of the phase flip channel and the energy exchanging channels. 
I.e, although there is information about the phase returning to the system via reversed 
jumps in the energy exchanging channels, the overall information flow about the phase 
is from the system to the environment due to the phase flips occurring in the dynamics. 
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This example illustrates, how in complex quantum dynamics it requires subtlety to 
determine whether memory effects play a role in the system dynamics. 

The NMQJ-approach gives a scheme to detect whether there are memory effects 
in the system dynamics although there are no revivals of coherence in the level of the 
density matrix. In [T7j the concept of non-Markovianity was clarified by defining a 
measure for non-Markovianity. It measures the amount of information flowing from the 
environment to the system. For the spin-boson model this measure is independent of the 
parameter e, i.e., the maximum value for non-Markovianity is obtained when there are 
no phase flips occurring in the dynamics. The NMQJ-description thus gives a similar 
picture of the memory effects in the dynamics of the spin-boson model as the measure 
of non-Markovianity. 

4. Conclusions 

We derived a time-convolutionless master equation describing the spin-boson model in 
the weak coupling limit. The equation describes the short time memory effects in the 
dynamics manifesting themselves through negative decay rates in the master equation. 

The memory effects are studied via the non-Markovian quantum jump approach. 
This description shows memory effects on the level of singe wave function realizations 
although there are no revivals of coherence in the density matrix description. This 
behavior demonstrates, how in complex quantum dynamics the co-operative action of 
different channels complicates the detection of memory effects in the dynamics. 
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